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The detailed investigation of electronic and magnetic properties of polymeric coordination ma- 
terials with accurate ab initio quantum mechanical methods is often computationally extremely 
demanding because of the large number of atoms in the unit cell. Moreover, usually the available 
structural data are insufficient or poorly determined especially when the structure contains hydro- 
gen atoms. In order to be able to perform controlled ab initio calculations on reliable structures, 
we use a two step approach to systematically prepare model structures for polymeric coordination 
compound systems and to relax them to their equilibrium configuration. First, a structure is con- 
structed on the basis of a crystallographic database and optimized by force field methods; in the 
second step, the structure is relaxed by ab initio quantum mechanical molecular dynamics. With 
this structure we perform accurate electronic structure calculations. We will apply this procedure 
to a Fe(II) triazole compound and to a coordination polymer of Cu(II) ions with 2,5-bis(pyrazol-l- 
yl)-l,4-dihydroxybenzene. 
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I. INTRODUCTION 

Polymeric complexes form a class of materials that at- 
tract interest from chemists and physicists. The pos- 
sibility to apply organic synthesis methods to produce 
families of compounds with similar structures but sub- 
tly differing physical behavior offers an opportunity to 
design materials with planned properties. 

Materials showing temperature or pressure induced 
high-spin - low-spin transitions are currently being dis- 
cussed for applications in data storage or spintronics. 
Such magnetic transitions can be optically triggered, 
raising hopes for technological applications. The co- 
operativity of such light-induced phase transitions also 
raises questions concerning the underlying microscopic 
processeai»2ijL4. Another important class of polymeric 
coordination compounds are, in general, low-dimensional 
magnetic materials where the coupling strengths between 
the metal ions can be tuned by chemically modifying the 
ligands and the transition metal ion coordination. Effects 
of, for instance, dimensional crossover can be studied in 
this way. 

Most of these compounds generally contain a large 
number of atoms in the unit cell and a very low sym- 
metry (in 80% of the cases triclinic or monoclinic) mak- 
ing the first principles study of the electronic and mag- 
netic properties of these systems under different condi- 
tions (pressure, temperature, light irradiation) computa- 
tionally very demanding. Moreover, since an accurate de- 
termination of hydrogen positions with X-ray diffractom- 
etry is very difficult, usually the available structure data 



are poorly determined implying unstable quantum me- 
chanical calculations. For some cases, the crystal struc- 
tures are not even known, due to absence of single crys- 
tals suitable for an X-ray structure analysis. In order 
to overcome these shortcomings, we use in this commu- 
nication a two-step approach which allows to perform 
accurate ab initio density functional theory (DFT) cal- 
culations on reliable model structures for coordination 
polymers. It should be noted, that the expression 11 ab 
initio prediction of crystal structures" is also used in the 
sense of" a priori prediction", i.e. predictions without ref- 
erence to diffraction data. Possible crystal structures of 
organic compounds with 20-60 atoms can be predicted 
within some days or weeks using force field methods 5 . 
This works also for organometallic compounds, e.g. pen- 
tamethylferrocene Fe(C5H5)[(C5CH3)5}&. Quantum me- 
chanical methods are much more time-consuming; hence 
their use is limited to local optimizations starting from a 
given structural model. Thus, quantum mechanical and 
force field methods can be combined by making first a 
global optimization using force field methods, followed 
by a local optimization with quantum mechanical meth- 
ods. This was done for many inorganic systems, e.g. 
NaCl or MgF^, and a few simple organic compounds, 
e.g. glycole C 2 H 4 (OH) 2 and glycerol C 3 H 5 (OH)^. In 
most other cases of crystal structure predictions, quan- 
tum mechanical calculations were only used for secondary 
calculations, e.g. to calculate molecular geometries and 
electrostatic charges, or to derive intramolecular or in- 
termolecular potential curves for subsequent force field 
optimizations^. 
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If experimental data are known, e.g. from crystal 
data of analogous structures, this knowledge should of 
course be used as far as possible in order to generate 
reliable model structures and to calculate solid state 
properties. Here, we apply ab initio DFT methods 
to understand the relationships between the structures 
and the electronic and magnetic properties of both i) 
spin crossover materials and ii) low-dimensional spin 
systems. We will consider two representatives from 
this large class of compounds: i) a Fe(II)triazole com- 
poun d 10 ' 11 which is thought to have a polymer structure 
and ii) a coordination polymer of Cu(II) ions with 2,5- 
bis(pyrazol-l-yl)-l,4-dihydroxybenzeneA2,, (we shall de- 
note it as [Cu(bpydhb)]„, alternative names given in the 
literature are CuCCB%. 

For the Fe(II)triazole syste m 10 ' 11 , there is no direct 
accurate determination of the crystal structure of the 
complex as it cannot be crystallized sufficiently well. 
A chain-like structure for such compounds of formula 
[Fe(Rtrz)3]A2.S where A denotes the anion was predicted 
from EXAFS spectroscopy^. For analogous Cu(II) com- 
pounds this chain structure was confirmed by single crys- 
tal X-ray analysis (see e.g. Refs. fl5lfl6l ). Polymeric iron 
triazole compounds always show a poor crystallinity, and 
single crystals cannot be grown. Hence the X-ray struc- 
ture analyses are limited to molecular compounds, e.g. 
trimers containing three Fe(II) ions connected by 6 tri- 
azole ligands and 6 terminal water molecules 17 . Poly- 
meric iron triazole complexes have been studied by pow- 
der diffraction; but apart from a preliminary indexing 
of [Fe(Htrz) 3 ](C10 4 ) 2 • 1-85 H 2 0±&, no structural data 
have been published so far. On the other hand, these 
polymeric compounds show interesting high-spin-low- 
spin transitions in the solid state (see e.g. Ref. Il9l). and 
quantum mechanical calculations could help to under- 
stand the mechanism of this spin crossover. Since all cal- 
culations of electronic and magnetic properties require 
a reliable and reasonable model of the crystal structure, 
we constructed the structures of a hypothetical polymer 
crystal. All available information of molecular iron tria- 
zole systems, and of analogous Cu(II) polymers was in- 
cluded. In order to make the ab initio DFT calculations 
feasible we constructed at first some models containing 
small triazoles (e.g. methyl-triazole) and simple anions 
(e.g. F~). Our rationale is that the model compound 
still has the same coordination and immediate environ- 
ment of the Fe(II) ion as the real material which should 
thus allow us to study the mechanisms of the high-spin 
- low-spin transition. Structures with more complex an- 
ions will be considered later. 

The second system that we study, the coordination 
polymer of Cu(II) ions with 2,5-bis(pyrazol-l-yl)-l,4- 
dihydroxybenzene, [Cu(bpydhb)]„, has been recently 
synthesized by Wagner et alM* The Cu 2+ ions in this 
compound form spin- 1/2 antiferromagnetic Heisenberg 
chains running along the z-axis. We will analyze the in- 
fluence on the electronic properties of slight structural 
modifications, like substitutions of small organic groups 



on the dihydroxybenzene or introduction of water ligands 
which change the coordination of Cu from square to oc- 
tahedral. The structure of Cu(II) with 2,5-bis(pyrazol- 
l-yl)-l,4-dihydroxybenzene has been obtained by X-ray 
powder diffraction^ and we will apply our two-step ap- 
proach for the modified structures based on these data. 



II. METHODS 

We derive and investigate model polymeric coordina- 
tion compound structures in two steps. First, we cre- 
ate a rough structure on the basis of crystallographic 
databases 2 ^. If no crystal structures of similar com- 
pounds are known - like in the case of Fe triazoles - 
we construct a hypothetical crystal structure having the 
highest possible symmetry, a small number of atoms per 
unit cell and a sensible arrangement of the individual 
fragments. The crystal structure is optimized by force 
field methods including the optimization of the unit cell 
parameters, but maintaining the crystallographic sym- 
metry. We used the Cerius 2 program package 2 ^ with a 
modified Dreiding 2 ^ force field; atomic charges were cal- 
culated by the Gasteiger method 2 -^. Afterwards we pro- 
ceed to improve the structure by relaxing it again, but 
now with an ab initio quantum mechanical molecular dy- 
namics method 2 ^. Note that we employ an all electron 
method already at this step in order to obtain the best 
possible equilibrium structure even for demanding transi- 
tion metal complexes. Finally, we analyze the properties 
of the relaxed model structure with an all electron ab 
initio DFT method of the highest precision 2 ^. 

Our Car Parrinello (CP) ab initio molecular dynam- 
ics (AIMD) calculations 2 ^ are performed with a projec- 
tor augmented wave (PAW) basis seiM. This is an all 
electron calculation based on the frozen core approx- 
imation. The exchange-correlation potential is evalu- 
ated in the generalized gradient approximation (GGA) 
as parametrized in Ref. \2E. After a first analysis of the 
electronic properties for the classically prepared model 
in order to determine the structural stability, we employ 
an AIMD with friction to refine the structure. We im- 
pose constraints on the atomic coordinates in order to 
preserve the space group symmetry of the compound. 

For the final analysis step we employ the full poten- 
tial linearized augmented plane wave (LAPW, APW+lo) 
method as implemented in the WIEN2k code 2 ^. 



III. RESULTS 

A. Fe(II)-Triazole 

The Fe(II) spin crossover molecular systems have been 
intensively studied in the past^S and various theories have 
been developed about the mechanism of the high-spin - 
low-spin transition. A few years ago, Kahn et alM pro- 
posed that the cooperativity of this transition may be en- 
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hanced by replacing molecular crystals by polymers, since 
in the latter the metallic sites are bridged by chemical 
linkers through which one may get efficient intersite cou- 
pling. A first principles DFT study of such systems would 
bring some light into the properties of these materials. 
Unfortunately, to date it has not been possible to get 
an accurate X-ray structure determination of Fe(II) spin 
crossover polymer systems, as discussed above, which is 
an indispensable input for ab initio DFT calculations and 
therefore, to our knowledge, no reliable ab initio studies 
were performed. With the two-step method, we have ob- 
tained - based on known data for molecular Fe(II) spin 
crossover systems - optimized model polymer structures 
which can be used as input for ab initio DFT calcula- 
tions. In Fig. Q] we present the structure of the poly- 
mer Fe(II) compound with 4-methyl-l,2,4-triazole moi- 
eties as bridging ligands between the Fe(II) cations. As 
anions we chose simple F~ ions in order not to compli- 
cate the DFT calculations. The structure has P2 1 /m 
symmetry; a higher symmetry was not possible with this 
chemical composition. As an initial structure with low 
spin character, the distances between the Fe(II) centers 
and the neighboring N atoms of the triazole rings were 
chosen to be e?Fe-N = 1-97 A, close to the Fe-N distance of 
c^Fe-N = 2.0 A that is known to be typical for the low spin 
variants of Fe triazole complexes 2 ^ while for the structure 
with high spin character we chose d^e-N = 

2.2 A . In Fig.[2] 

we present the spin-polarized density of states (DOS) ob- 
tained 26 - for this system after structural optimization for 
the low spin (Fig. [2] (a)) and the high spin (Fig. [2] (b)) 
states. Shown is the total Fe d DOS for majority (upper 
panel) and minority (lower panel) spin states. Since Fe 
is in an octahedral environment of N, the 3 d states split 
into lower t 2g (d xy ,d xz ,d yz ) and higher e g (d x 2_ y 2 : d z 2) 
states. In Fig. [2] (a) the 6 electrons of Fe(II) 3d 6 occupy 
fully the t2 g states (see DOS manifold between -2 and 
-0.2 eV) for both spins thus corresponding to a low spin 
S=0 state. The e g states are completely empty and lie 
above the Fermi level. In Fig. [2](b), the majority (spin 
up) t2 g and e g states are occupied (upper panel DOS in 
the range between -3.8eV and -0.6eV) while the minority 
(spin down, lower panel) t 2g states are about one-third 
occupied and two-thirds empty. The minority e g states 
are completely empty (states in the range of energies be- 
tween leV and 3eV, lower panel). These results indicate 
that the system is in a high spin S=2 state. A detailed 
investigation of the magnetic properties of this polymer 
system for both spin states will be presented elsewhere 2 ^ 



B. Cu(II)-Polymer 

The Cu(II) coordination polymer [Cu(bpydhb)]„ syn- 
thesized by Wagner et al^ crystallizes in the triclinic 
space group PI (No. 2) and consists of stacks of chains of 
deprotonated 2 ,5-bis (pyrazol- 1-yl)- 1 ,4-dihydroxybenzene 
molecules connected by Cu 2+ ions. This system behaves 
as a homogeneous antiferromagnetic spin-1/2 Heisenberg 




Figure 1: Two unit cells of the triazole model polymer. The 
usually complicated anions have been replaced by single atom 
F~ anions. 



chain with an exchange coupling constant J = 21.5 K/ke 
estimated from magnetic susceptibility measurements^. 
The modular nature of this polymer makes it an adequate 
candidate to study chemically feasible structural modi- 
fications which may change the electronic and magnetic 
properties. Of special interest is the effect of dimension- 
ality, i.e. tuning the crossover from one dimension to 
higher dimensions with the inclusion of adequately mod- 
ified ligands. Since in this communication we want to 
stress the procedure we follow to design these structures 
and how we calculate the electronic properties, we will 
concentrate here on one substitution, namely two hydro- 
gen atoms in the hydroquinone ring by two amino groups 
and show the efficiency of the two-step process. A sec- 
ond issue in this context that requires some attention is 
the study of the effect of two additional ligands at the 
Cu(II) ion, since such an octahedral coordination of cop- 
per is also frequently observed. We will here present the 
relaxation process on [Cu(bpydhb)]„ with two additional 
water molecules in the environment of copper. 

In Fig. [3] we present the relaxed structure of 
[Cu(bpydhb)]„ where two hydrogen atoms in the hy- 
droquinone rings have been substituted by two amino 
groups. This structure was obtained first by applying the 
force field method, which already has the correct informa- 
tion of the chemical distances and angles from data bases. 
A quantum mechanical evaluation of the forces on the 
atoms (see Table [J) for this structure, performed with the 
full potential LAPW method 2 ^ (FPLAPW) still shows 
an unstable structure. Our aim is to relax this structure 
showing forces in the hundreds of mRyd/ae with quan- 
tum mechanical molecular dynamics calculations, AIMD, 
until we find forces smaller than 25 mRyd/ae which we 
set as upper bound for the force values. In order to get 
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Figure 2: Spin polarized Fe d density of states 
of the Fe triazole structure Fig. [Tj in units of 
states/electronvolt/atom/spin. (a) for the S=0 state of Fe 
(b) for the S=2 state of Fe. The upper and lower panels of 
the figure show respectively the density of states for the ma- 
jority and minority spin states. The Fermi energy is indicated 
by a dotted line. 




Figure 3: Four unit cells of [Cu(bpydhb)] n (see text for the 
detailed formula) where two hydrogen atoms in the hydro- 
quinone rings have been substituted by two amino groups. 
We find that the relaxed structure shows some tilting of the 
amino groups out of the plane defined by the benzene ring. 



Table I: Forces^ between the atoms for the Cu(II)-NH2 poly- 
mer structure designed with the force field method. Note that 
the forces were determined under the constraint of inversion 
symmetry; thus there are no forces acting on the Cu(II) ion 
center. Units are in mRyd/aB. See Table lHII for identification 
of the atom location. 



Atom 


Fx 




F z 


Cu 














-47.4 


-62.7 


51.3 


Nl 


89.1 


32.8 


-42.6 


N2 


22.4 


27.2 


-95.8 


N3 


-94.2 


-66.8 


33.5 


CI 


35.0 


22.9 


1.8 


C2 


-92.7 


-29.2 


111.1 


03 


-39.2 


-41.7 


294.8 


C4 


-239.7 


-104.1 


16.5 


Co 


-26.6 


-12.8 


-304.2 


C6 


106.7 


88.9 


46.0 


HI 


128.6 


-28.7 


215.0 


H2 


255.3 


106.1 


-34.2 


H3 


23.8 


25.6 


-279.0 


H4 


4.8 


10.8 


-18.3 


H5 


-16.2 


64.5 


56.3 



structures with forces smaller than 1 mRyd/aB within the 
FPLAPW method, a fine tuning between the various first 
principles DFT methods we employ has to be considered. 
For our purposes here, a benchmark of 25 mRyd/aB is 
reasonable. 

In Table HI] the forces after the AIMD relaxation are 
shown and we observe that the structure is relaxed to 
our required precision. The fractional coordinates of this 
structure are given in Table HTT1 where the symmetry and 
cell parameters were kept unchanged with respect to the 
unsubstituted Cu polymer system [Cu(bpydhb)]„— a = 
5.1723 A, b = 7.9587 A, c = 8.2298 A, and angles a = 
118.221°, (3 = 91.520°, 7 = 100.148°. 

In a different modification of the Cu(II) polymer struc- 
ture, we want to investigate the influence of a change of 
the Cu coordination from four to six on the electronic 
properties of the Cu(II) polymer. We added two water 
molecules bound to the Cu center and proceeded as out- 
lined for the case of the NH2 substitution: A first struc- 
ture was obtained from the force field method, yielding 
high forces when we controlled it with an FPLAPW cal- 
culation. Further relaxation with AIMD yields forces 
around 20 mRyd/as- The corresponding structure is 
shown in Fig. [4j and the fractional coordinates are given 
in Table IIVI The new unit cell is triclinic with space 
group PI (No. 2) as in the original Cu(II) polymer. 
The cell parameters are a = 5.2341 A, b = 11.2493 A, 
c = 8.0721 A and the angles are a = 117.6108°, = 
68.8218°, 7 = 127.1551°. Note the variation of these pa- 
rameters with respect to those of [Cu(bpydhb)],r^, since 
in the relaxation process of the new structure the unit cell 
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Table II: Forces^ between the atoms for the relaxed Cu(II)- 
NH2 polymer. Units are mRyd/aB. 



Atom 


Fx 




F z 


Cu 














12.7 


16.9 


21.7 


Nl 


3.9 


0.7 


-17.2 


N2 


12.5 


-2.3 


8.4 


N3 


-21.0 


-9.4 


-5.5 


CI 


-1.0 


-9.3 


-19.4 


C2 


10.3 


3.6 


-2.1 


C3 


-7.1 


-7.0 


14.6 


C4 


-10.6 


-12.7 


-3.1 


Co 


-12.0 


-6.0 


-19.5 


C6 


12.0 


3.4 


-0.7 


HI 


1.4 


-0.7 


3.5 


H2 


3.1 


-1.7 


-1.4 


H3 


-0.4 


-0.3 


-4.1 


H4 


0.08 


-5.5 


-0.4 


H5 


1.1 


-0.9 


0.5 



Table III: Fractional atomic positions of nonequivalent atoms 
in Cu(II)-NH2 obtained after relaxation. The corresponding 
unit cell parameters are given in the text. 



Atom 


X 


y 


z 


Cu 


0.5 


0.5 


0.5 


O 


0.45952338 


0.34384463 


0.62047714 


Nl 


0.87129861 


0.67649196 


0.85703667 


N2 


0.81592808 


0.67403219 


0.69273582 


N3 


0.89342131 


0.81953377 


0.25313235 


CI 


0.48379900 


0.42607948 


0.80332348 


C2 


0.68236558 


0.58873536 


0.92870320 


C3 


0.03998649 


0.76736537 


0.66677421 


C4 


0.23751744 


0.83563448 


0.81694695 


Co 


0.12514376 


0.77281287 


0.93447831 


C6 


0.70113966 


0.66249861 


0.12137891 


HI 


0.20669437 


0.78002257 


0.06146726 


H2 


0.43843580 


0.91660682 


0.83522646 


H3 


0.04936408 


0.77481595 


0.53807198 


H4 


0.82019516 


0.86915420 


0.37967013 


H5 


0.94873662 


0.93065156 


0.22239484 



Table IV: Fractional atomic positions of nonequivalent atoms 
in Cu(II)-H20 obtained after relaxation. The corresponding 
unit cell parameters are given in the text. 
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Nl 


0.88662263 


0.72926141 


0.81944391 


N2 


0.82847386 


0.71119262 


0.65595860 


HI 


0.23644016 


0.90962503 


0.01364916 


H2 


0.47761970 


0.07047826 


0.76570604 


H3 


0.05611367 


0.85699006 


0.49132676 


H4 


0.20991897 


0.21421176 


0.81558567 


Ho 


0.02924495 


0.44665601 


0.79637256 


H6 


0.92725690 


0.45505960 


0.64249638 




Figure 4: One unit cell of the [Cu(bpydhb)] n polymer with 
water molecules bonding to the Cu(II) ion. The unit cell 
has been relaxed by the force field method but kept fixed in 
the AIMD relaxation. It differs considerably from the original 
[Cu(bpydhb)]„ polymer unit cell. The equilibrium positions of 
the water molecule lead to a nearly octahedrally coordinated 
Cu(II) ion. 



was allowed to relax properly in order to accommodate 
the water molecules. 

We investigate now the electronic properties of the two 
model Cu(II) polymers with a precise FPLAPW calcula- 
tion. A detailed account of the electronic and magnetic 
properties of both systems is presented elsewhere^. Here 
we would like to stress the effect of a planar Cu environ- 
ment (for the Cu-NH 2 case) in comparison of an octa- 
hedral environment (for the Cu-H 2 case) on the Cu d 
density of states. In Figs. [5] (a) and (b) the Cu d density 



of states for both modifications is presented. While for 
both cases Cu is in an oxidation 2+ with the d x 2_ y 2 or- 
bital half filled, the effect of the octahedral crystal field 
shows a splitting into ti g (d xy ,d xz ,d yz ) for energies be- 
low -1.5 eV and e g (d z 2,d x 2_ y 2) states above -1.5 eV (see 
Fig. [5] (b)) while in the planar environment all d states 
are split and the states contributing to the bonding at 
the Fermi surface are mainly of d x 2_ y 2 character in the 
local coordinate frame of Cu which is defined with the 
local z direction pointing from the Cu to the O of the 
H2O molecule and the y direction pointing from the Cu 
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dimensional polymer [Cu(bpydhb)]„ with small struc- 
tural modifications and a model Fe(II) triazole polymer 
using classical and ab initio DFT methods. We demon- 
strate that the use of a stepwise approach to the study of 
complicated coordination polymer materials is an effec- 
tive procedure to obtain reliable structures for accurate 
quantum mechanical analysis. We note that an exclusive 
use of ab initio molecular dynamics is very time con- 
suming while the combination of the classical with the 
quantum approaches speeds up the relaxation process 
considerably. We first employ a classical force field to 
create and relax plausible initial structures. In a second 
step, we proceed with ab initio molecular dynamics to re- 
lax these structures to a stable configuration. Finally we 
extract the electronic properties of the transition metal 
compounds with precise FPLAPW calculations. 



Figure 5: (a) Cu d density of states of the Cu polymer 
[Cu(bpydhb)]n with two NH2 groups, and (b) of the Cu poly- 
mer [Cu(bpydhb)] n with two water ligands. 



to the in-plane atom. 



IV. CONCLUSIONS 

We have investigated the structural properties and 
its effects on the electronic properties of the quasi-one- 
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